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Abstract 

iBeniamini. Burdzv. and Chenl ( 200?! ) introduced the notion of a shy coupling: a couphng of a Markov 
process such that, for suitable starting points, there is a positive chance of the two component processes 
of the coupling staying at least a given positive distance away from each other for all time. Among other 
results, they showed that no shy couplings could exist for reflected Brownian motions in bounded 
convex planar domain s whose boundar i es co ntain no line segments. Here we use potential-theoretic 
methods to extend this lBeniamini et al.l ( 200?! ) result (a) to all bounded convex domains (whether planar 
and smooth or not) whose boundaries contain no line segments, (b) to all bounded convex planar domains 
regardless of further conditions on the boundary. 

MSC 2000 subject classification: 60J65 



1 Introduction 



Motivated by the use of reflect ion couplings for reflected Brownian motion in a number of contexts (for 
ex ample efficient coupling as in Burdzv and Kendal^ 200cl, and work related t o the hotspot conjecture, a s 
in]Atar and Burdzv, '2OO2', and recent work on Laugesen-Morpurgo conjectures, Pascu and Gageonea , 20081 ), 
[Beniamini. Burdzv, and Chen (200^) introduced the notion of shy coupling and studied it in the contexts of 
Brownian motion on graphs and reflected Brownian motion particularly in convex domains. Shy coupling 
of two random processes occurs when, for suitable starting points, there is a positive chance of the two 
component processes of the coupling staying at least a given positive distance away from each other: this of 
course is in direct contrast to the more usual use of couplings, in which the objective is to arrange for the 
two processes to meet. I nterest in shy couplings is therefore focused on characterizing situations in which shy 
coupling cannot occur. IBeniamini et al.l (|2007l ) show non-existence of shy coupling for reflected Brownian 
motions in bounded convex planar domains whose boundaries contain n o line segments . The present note 
(influenced by techniques used to study non-confluence of F-martingales in lKendal] , ll99Cll ) uses rather direct 
potential-theoretic methods to obtain signiflcantly stro nger results in the case of reflected Brownian motion 
in convex domains (Theorems |7| and |8| below) . As with iBeniamini et al.l (|2007l ), the principal purpose is to 
contribute to a better understanding of probabilistic coupling. 

To flx notation, we begin with a formal discussion of coupling and shy coupling. In the following D is 
a measurable space, and the distribution of a Markov process on D is specifled as a semigroup of transition 
probability measures {Pf : t ^ 0,z e D}. (The semigroup-based deflnition allows us to take account of 
varying starting points, a typical feature of coupling arguments.) 
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Definition 1 (Co-adapted coupling). A coupling of a Markov process {P^ : t ^ 0,z G D} on a measurable 
space D is a family of random processes (X,Y) on D^, one process (X,Y) for each pair of starting points 
(^'OiDo) G D^, such that X and Y each separately are Markov but share the same semigroup of transition 
probability measures {P^ : t ^ 0,z e D}. Thus for each s,t ^ and z G D and each measurable A C D we 
have 

P [Xs+t G A I Xs = z, Xu : u ^ s] = P^{A) , 
P [Ys+t G A I Ys = z, Yu : s$ u s] = Pt"(A) . 

The coupling is said to be a co-adapted coupling if the conditioning in ([T]) can in each case be extended to 
include the pasts of both X and Y: for each s, t ^ and z G D and each measurable A C D 

P[X3+te AIX3 =z, Xu,Y^:0^u^s] = Pt"{A), 
P[Ys+te AIYs^z, Xu,Y^:0^u^s] = P^^IA) . 



In contrast to iBeniamini et alJ (|2007l )'s notion of Markovian coupling, we do not req uire (X, Y) to be 



Marko v. (This generalization is convenient but unimportant.) Note that the couplings in IBeniamini et al 
are all Markovian and hence co-adapted. 
We say that (X,Y) begun at (xo,yo] couples successfully on the event 

[Xt — Yt for all sufficiently large t] . 

Note that a co-adapted coupling can be adjusted on the simpler event [Xt = Yt for some t] so as to couple 
successfully on the new event; this need not be the case for more general couplings. In the remainder of this 
no te the term "coup l ing" w ill always be short for "co-adapted coupling". 



Beniamini et al.l (|2007l )'s notion of shy coupling for a Markov process on D is primarily concerned with 



the cases of Brownian motion on graphs and reflected Brownian motion on domains in Euclidean space. 
However their definition is expressed in general terms: suppose that D is actually a metric space equipped 
with a distance dist which furnishes a Borel measurability structure: 

Definition 2 (Shy coupling). A coupling {X,Y) is shy if there exist two distinct starting points Xg ^ yo 
such that for some £ > 

P[dist(Xt,Yt) > £ for all t 1 Xq =Xo,Yo ^-yo] > 0. (3) 

In words, X, Y has positive chance of failing to £-couple for some £ > and some pairs of starting points. 
We say that (X,Y) begun at {xo,yo) e-couples (for £ > 0) on the event 

[dist(Xt, Yt) — £ for some t] . 

We focus on reflecting Brownian motion in a bounded convex domain in finite-dimensional Euclidean space: 

Definition 3 (Reflecting Brownian motion). A reflecting Brownian motion in the closure D of a bounded 
convex domain D C M^, begun at Xq G D, is a Markov process X in D solving 



ft 

Xt = Xq + Bt + 



^{Xs)dL^ (4) 



Jo 

where B is n-dimensional standard Brownian motion, L''^ measures local time of X accumulated at the 
boundary 9D, and -v is a choice of an inward-pointing unit normal vectorfleld on 9D. 
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Unique solutions of Q exist in the case of bounded convex domains ( Tanakal . 19791 , Theorem 3.1); this is a 
consequence of results about the deterministic Skorokhod's equation. Note that reflected Brownian motion 
can be defined as a semimartingale for domains which are muc h more general than convex domains (see for 
example the treatment of Lipschitz domains in iBass and Hsu , 19901) ; however the results of this paper are 
concerned entirely with the convex case. 

Recall that a convex domain in Euclidean spaces can be viewed as the intersection of countably many half- 
spaces. Consequently the inward-pointing unit normal vectorfield "v is unique up to a subset of 9D = D \ D of 
zero Hausdorff (n— l)-measure. The local time term L''^ can be defined using the Skorokhod construction: it 
is the minimal non-decreasing process such that the solution of (|4]) stays confined within D. It is determined 
by the choice of normal vectorfield "v. 

Arguments from stochastic calculus show that all co-adapted couplings of such reflected Brownian motions 
can be represented in terms of co-adapted couplings (A, B) of n-dimensional standard Brownian motions, 
which in turn must satisfy 



At = 



dBs + 



dC. 



where C is an n-dimensional standard Brownian motion independent of B, and JJ, K are (n x n) -matrix- valued 
random processes, adapted to the filtration generated by B and C, and satisfying the identity 



+ ] 



(the identity matrix on 



A proof of this fact can be found in passing on page 297 of lEmervl (|2005l ). An explicit statement and a 
slightly more direct proof is to be found below at Lemma [6l 
Thus we will study X and Y such that 



Xt = xo 



Yt = 

JJ7 Jt + ] 



Bt + 





^{X,)dL 



X 
s J 





dB, 



[T d Cs + 



^(Ys)dq 



(5) 



The results of this note concern all possible co-adapted couplings, of which two important examples are: 

Ytll, so that dA is the reflection of dB in 



1. the reflection coupling: 

K = and Jt = l'"^' -2ete^r, where et = (Xt - Yt)/|lXt 
the hyperplane bisecting the segment from Xt to Yt; 



and its opposite, 

2. the perverse coupling: 

K = and Jt = +2 et et'^, so that the distance ||X — Y|| has tra j ectori es purely of locally bounded 

variation. (This is related to the uncoupled construction in Emery, 1989t Exercise 5.43 of the planar 
process Z" with deterministic radial part.) 

In practice, as is typically the case when studying Brownian couplings, general proofs follow easily from the 
special case when K = an d J is an orthogonal matrix. (Heuristic remarks related to this observation are to 
be found in lKendalll , I2OO7I , Section 2.) 
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Note that the existence question for reflection couphngs in convex domains is non-trivial (|Atar and Burdzv 



2004] estabhsh existence fo r the more general case of lip domains), but is not relevant for our purposes. 



Beniamini et alJ (j2007l ) use ingenious arguments based on ideas from differential games to show that 



a bounded convex planar domain cannot supp ort any shy couplings of reflected Brownian motions if the 



boundary is and contains no line segments (jBeniamini et al 



20071 , Theorem 4.3). Here we describe the 



use of rather direct potential-theoret ic methods (summarized in section [2]); in section [3] these are used to 
generalize the iBeniamini et al.l ([2007j) result to the cases of (a) all bounded convex domains in Euclidean 
space of whatever dimension with boundaries which need not be smooth but must contain no line segments 
(Theorem [7|), and (b) all bounded convex planar domains (Theorem [8}. Further extensions and conjectures 
are discussed in section [H 



2 Some lemmas from stochastic calculus and potential theory 

The following two lemmas from probabilistic potential theory are fundamental for our approach. 

Lemma 4. For D C a bounded domain, for fixed e > 0, consider the closed and bounded (and 
therefore compact) subset o/M" x given by 

F = (D X D)\{(x,i)) :dist(x,y) < £}. 

Suppose there is a continuous function W : ¥ —>■ such that the following random process Z is a 
supermartingale for any co-adapted coupling of reflecting Brownian motions X and Y in the closure 
D, when S is the exit time of (X,Y) from F; 

Zt = ¥(XtAs,YtAs)+tAS. 

Then for any such coupling almost surely S < oo and dist(Xs, Ys) — e. 

Proof. Being a continuous function on a compact set F, W must be bounded. Thus for any co-adapted 
coupling (X,Y) the supermartingale Z is bounded below, and so almost surely it must converge to a flnite 
value at time oo. But boundedness ofW also means that Z must almost surely tend to inflnity on the event 
[S = oo]. These two requirements on Z force the conclusion that P [S < oo] = 1. Since X and Y each reflect off 
9D, it follows that the exit point (Xs, Ys) must belong to the part of 9F which is contained in the boundary 
of {(x,^) : dist(x,i|) < e}; hence almost surely dist(Xs,Ys) = £. □ 

Consequently, if one can exhibit such a W for a specified D then any coupled X, Y in D must £-couple. Indeed 
the existence of s uch a W implies a unifor m bound on exponential moments of the £-coupling time S. 



As implied by Beniamini et al. ( 200?!, Example 4.2), the obvious possibility V = c||X — Y|| " (for positive 



c, a with a small) is not suitable here; the perverse coupling of the previous section is an example for which 
all such "'i' lead to ^(XtAS , YtAS ) + t A S being a strictly increasing process when neither X nor Y belong 
to 9D. Nevertheless we will see in the next section how to construct W which work simultaneously for all 
co-adapted couplings in a wide range of bounded convex domains. 



Before embarking on this, it is convenient to introduce a further lemma motivated bv lltd (|l975l )'s approach 
to stochastic calculus using modules of stochastic differentials: if Z = M-l-A is the Doob-Meyer decomposition 
of a semimartingale Z into local martingale M and process A of locally bounded variation, then write 
Drift dZ = dA and write (dZ)^ — d[M, M] for the differential of the (increasing) bracket process of M. 
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Lemma 5. The conclusion of Lemma^ holds if there is a continuous O : F ^ M such that, for all 
co-adapted couplings (X,Y) as above, 

Z = <1)(X,Y) is a semimartingale; 

moreover the stochastic differential dZ satisfies the following random measure inequalities (with 

a,h> 0): 

1. (Volatility bounded from below) (dZ)^ > adt up to time S; 

2. ( Drift bounded from above ) Drift d Z < b d t up to time S . 

Proof. For c, A > to be chosen at the end of the proof, set ¥(x,y) = c(l — exp(— A<D(x,y))) and apply Ito's 
lemma to V(X, Y) = c(l — exp(— AZ]] up to the exit time S: 

Drift d¥(X,Y) = cAexp(-AZ) (Drift dZ - |A{dZ)2) <^ -f A exp(-AZ) (Aa - 2b) dt 

(where the inequality is viewed as an inequality for random measures). Fix A > 2b/Q and c > 2A~ ^ao" 2b '^'^ 
The above calculation then shows that ^(XtAS,YtAs) + (t AS) is a supermartingale, so the conclusion of 
Lemma [H holds. □ 

As mentioned above in section [ Tl the following stochastic calculus result is to be found unannounced and 



m passmg on page 297 of lEmervl POOSI ). It provides an explicit representation for co-adaptively coupled 



Brownian motions. We state it here as a lemma and indicate a direct proof in order to establish an explicit 
statement of the result in the literature. 

Lemma 6. Suppose that A and B are two co-adapted n- dimensional Brownian motions. Augmenting 
the filtration if necessary by adding an independent adapted n- dimensional Brownian motion D, it is 
possible to construct an adapted n- dimensional Brownian motion C, independent o/B, such that 



A = 



dB + 



dC, 



where JJ, K are (n x n) -matrix-valued predictable random processes, satisfying the identity 

J + IK^ K = l'^' (the identity matrix on R"^) . 

Proof. Certainly the quadratic covariation between the vector semimartingales A and B may be expressed 
as dA dB^ = JJ^dt for a predictable (n x n) -matrix- valued predictable random process JJ such that the 
symmetric matrix inequality ^ JJ^ J ^ l'"^' holds in the spectral sense (0 ^ Jx ^ x^x for all vectors 

x); this is a consequence of the Kunita-Watanabe inequality. Since JJ^ J is a symmetric contraction matrix 
we find (I^ J)'' ^ Hi as k — > oo, where Hi is the orthogonal projection onto the null-space of JJ — l'"^'. 
This exhibits Hi as a measurable function of J. We may now extract Aj — sup{x^(J^ J— Hi)x : ||x|| — 1} as 
a measurable function of J, renormalize and study (J^ JJ— Hi)^/A2'^ H2, and continue so as to represent 

N 

JJ^J = Hi-F^A^Hi, (6) 

1=2 
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where Hi, H2, . . . , Hn are disjoint orthogonal projections, N ^ n, and 1 > A2 > . . . > An > and all 
quantities are measurable functions of J. 

Set Ho = l'"^' -Hi - H2 - . . . - Hn and define 



K = Ho + ^^l-A^Hi 
1=2 

as the non-negative symmetric square root of l'"^' — J. The above spectral approach defines K as a 
measurable function of the matrix J, and hence we may view K as well as J as predictable (n x n)-matrix- 
valued predictable random processes, now satisfying JJ^ J + K^ K — l'^'. 

Moreover we may construct a pseudo-inverse of K as a further predictable symmetric process: 



i=2 



so that K+ K = K K+ = l'"' -Hi. 
Now define C by 

dC = K+(dA- J^dB) +HidD. 
The quadratic variation of d A — d B is 

(dA- jTdB){dA- J^dB)^ = (I'^^'-J^Ddt^K^Edt, 

and so the stochastic differential K^(d A — JJ^ dB) has finite quadratic variation (l'"^' —Hi) dt and therefore 
(by the theory of stochastic differentials of continuous semimartingales) it is a martingale differential. 
Hence C is a continuous martingale with quadratic variation d C d = I*"^' d t, by which we may deduce that 
C is n-dimensional Brownian motion. Moreover dCdBT=K+(dA-J^dB)dB^=K+(J^-J^)dt = 0, 
so C is independent of B. 

Finally K^dC = (l'^' -Hi](d A - dB) (use Hi = 0) and Hi{dA- J^dB) is a martingale differ- 
ential with quadratic variation 

Hi(dA - J^dB)(dA - jTdB)^Hi = Hi(l'^' - J[)Hi dt = HiK^KHidt = 0. 
So d C = d A - d B, which establishes the result. □ 



Emervl (|2005l ) uses an approximation argument, which we circumvent here by exhibiting IK as a predictable 



process via the spectral decomposition ([6}. Note also that the extra Brownian motion D is not required if 

JJ^ J < l'"^' , in which case Hi = 0. 



3 Theorems and proofs 



The first theorem generalizes the planar result of iBeniamini et al.l (|2007l , Theorem 4.3) to the case of higher 



dimensions, and removes the requirement of boundary smoothness. 

Theorem 7. Let D be the closure of a bounded convex domain in such that 9D contains no line 
segments. Then no co-adapted coupling of reflected Brownian motions in D can be shy. 
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Proof. It suffices to exhibit, for any fixed £, a function O satisfying the two requirements of L emma HI 
Moti vated by a similar construction used to establish the convex geometry of small hemispheres (KendalJ, 
199ll ), define (for p ^ D and 6 > 0) 



Vp,6(x,y) 



Illy 



(7) 



Thus Vp_6 is a hyperbolic perturbation of |||x — y|p based on the pole p. For any fixed e > and for all 
sufficiently small 6 > depending on £, p, and the geometry of D we show that O — Vp ^ satisfies the 
requirements of Lemma [5l The result then follows by applying Lemmas [5] then HI 

The first step is to establish Lemma [5] requirement [1] (volatility bounded from below). Fixing £ > 0, 
suppose by virtue of Lemma [6] that X, Y satisfy ([5| for some co-adapted matrix processes JJ, K. Applying 
Ito's lemma and the fact that B and C are independent n-dimensional standard Brownian motions, 

(d(D{X,Y))2 = (||X-Y + 6(X-p) -J(X-Y + 6(Y-p))||2 + ||]K(X-Y + 6(Y-p))||2)dt 

^ ||X-Y+6(X-p)-J(X-Y + 5(Y-p))||2dt 

^ (||X-Y+5(X-p)||-||X-Y+6(Y-p)||)2dt, (8) 

since the third equation of ([5| implies that the linear map J is a contraction. Suppose now that ||X — Y|| > £. 
If (for example) 



5 < 



2 sup{dist{p,w) : w e D} 



(9) 



then 



IX- Y+ 6fX 



= ||X-Y+6(Y-p] +6(X-Y)|12 

= ||X-Y+5(Y-p)||2 + 25||X-Y||2 + 252(X- Y, ^ - p) 

P) 

P) 



|X- Y + 6(Y-p)||2 + 26||X- Y||2 -252|(X- Y,^ 



= ||X-Y + 6(Y-p)||^ + 26||X-Y|lM 1-5 



X-Y X+Y 
||X-Y||' 2 



IIX-YII 
IIX- Y- 



6(Y-p)||2 + 6£2. (10) 



This in turn implies 



|X-Y+6(X-p)||- 



X-Y+6(X-Y)|| ^ 



6£2 



lX-Y+5(X-p)|| + ||X-Y+6(X-Y)|| 



(2 + 6) diamD + 6 sup{dist(p, w) : w e D} 



(11) 



Thus for all small enough 6 > (bounded by ([9])) there is a constant a — a(£,6,p,D) > such that 
d(I)(X,Y))^ > adt while ||X — Y|| > £, no matter what co-adapted coupling is employed. This establishes 
requirement [l] of Lemma [H Note that we have not yet used the condition that 9D be free of line segments. 
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The second step is to establish Lemma [5] requirement [2] (drift bounded from above). Prom fSj), 



Drift d(D(X,Y) = 

(X - Y + 5(X - p),t(X)) dL^ + (Y - X - 6(Y - p),t(Y)) dL^ + (n - | tr(JI + J^) dt . (12) 

Since tr(J + J^) ^ — 2rL, requirement [2] of Lemma [5] is estabhshed with b = 2n if we can show, for x, ij in D 
with j|x — 1)11 ^ £, 

(x — y + 6(x — pjj'vfx)) < when X <E 8D , 
(y -X- 6(ij -p),-^(y)) < whenye9D. 

Both inequalities hold for all small enough 6 > depending on e, p, and the geometry of D. Consider for 
example the second inequality in (fTSj) . The Euclidean set 

K — {(x,y,v) : x e D,y e 9D,v is unit inward-pointing normal to 9D at y}, 

is closed and bounded hence compact. By convexity of D the inner product (y — x,v) is non-positive on the 
subset K and can vanish on K only when x = y or the segment between x and y lies in 9D. Moreover this 
inner product is a continuous function on K. But 9D contains no line segments and so (y — x, v) is negative 
everywhere on the compact set K\{(x,y, v) : ||x— y || < e}, and therefore must satisfy a negative upper bound. 
Thus the second inequality of (fTSjl also holds when ||x — y || ^ e, for all small enough 6 depending on e, p, 
and the geometry of D . The first inequality follows similarly. 

Hence both requirements of Lemma [5] apply for O = Vp g for any fixed e > once 6 is small enough 
(depending on e and the geometry of D). It follows from Lemma H] that any co-adapted coupling X,Y 
eventually attains dist(X,Y) ^ e for any £ > 0, and hence no co-adapted coupling can be shy. □ 

Theorem 8. Let D foe a bounded convex planar domain. Then no co- adapted coupling of reflected 
Brownian motions in D can be shy. 

Proof. The proof strategy is the same as for Theorem [7l exhibit, for any fixed e > 0, a function O satisfying 
the requirements of Lemma [5l However the function will now be the minimum CD = Oi A . . . A Ok of 
a finite number of fu nctions Ot of the form Vp_6 or a mild generalization thereof. The Tanaka formula 



([Revuz and Yorl , ll99ll , VI Section 1.2) shows that if each satisfies requirement [T] of Lemma[5l and satisfies 
requirement [2] on Oi = O, then O satisfies both requirements and so we may apply Lemmas [5] and then [4] as 
above. 

The mild generalization modifies Vp_5(x,y) by a small perturbation. This is required in order to deal 
with the possibility that 9D contains parallel line segments. 

Vp,6.s(x,y) = i||x-y||2 + f (||x-p||2-||y-p||2) 

= i||x-y||2 + K(x-y,2S±V-p), 

where now K = 6exp (-|||^^ -p||) . (14) 

Thus the asymmetry of the hyperbolic perturbation Vp_5_s now depends on the distance of the mid-point 
from the pole p. 
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The first step is to establisli Lemma [5] requirement [T] (volatility bounded from below). The argument 
runs much as for Theorem [7l but is further complicated by the need to deal with parallel line segments in 
9D. For convenience we introduce the notation 

Z = ^-P, e = Z/llZll. (15) 

Thus (using the fact that J is a contraction) 

(dVp,6,s{X,Y))Vdt = 

||X-Y+k(X-p-^(X-Y,Z) e)-J(X-Y+K(Y-p + ^(X-Y, Z)e))\\^ + 
+ II K(X - Y + k(Y - p + ^(X - Y, Z> e))||2 

^ (||X-Y+K(X-p-^(X-Y,Z)e)||-||X-Y+K(Y-p+^(X-Y,Z)e)||)'. 
Now consider the difference of the squared norms of the summands: if 

S > sup{dist(p,w) : w e D}, (16) 
K < 5 < e/sup{dist(p,w) : w e D} (17) 

(so that in particular 1 — > 0) then 

||X-Y+K(X-p-^(X-Y,Z)e)||2-||X-Y+K(Y-p + ^(X-Y,Z)e)||2 

= 2k(X- Y+ kZ,X- Y- |(X- Y,Z) e) 

= 2k (||X- Y||2 - m(x- Y,e)2 + k||Z|| (1 - ^)(X - Y, e> 

^ 2k||X-Y||2(i-^) _2K(K||Z||(1-^)|(X-Y,e) 

>2K||X-Y||^(l-M) (l-K||Z||K^) 

^ 2i:||X ^"Ptdi^t(P;>v)^>veD} ^ ^1 ^, sup{dist(p,w) :wg D} ^ ^ 

It follows from inequalities p6|) and (fTTjl that we obtain a positive lower bound on 

||X-Y+K(X-p-^(X-Y,Z)e)||-||X-Y+K(Y-p + ^(X-Y,Z)e)|| 

subject to the further condition (required to obtain a lower bound on each of the two norms in the difference 
above) 

K s; 6 < e/ (sup{dist(p,w) : w e D}+ |diam(D)) . (18) 
Now, for a finite set of poles p^, to be determined below in the "localization argument", consider 

TR 

(S>[x,v) = A (^pr.s.sKy) Av^+_5_s{y,'') Av^-_5_s(^,^j) Av^-_6_s(y,'')) • (i9) 

i=l 

Subject to S and k satisfying the conditions ([l6l [iTl [18} as p runs through the poles it then follows from 
the Tanaka formula that O satisfies requirement [T] of Lemma [H 
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Figure 1: Circle C centred in D of large radius R, with intersection points produced by extending the maximal 
line segments in 9D (for simplicity D is chosen to be a polygon). Note that £i and £3 are parallel. 

The second step is a "localization" argument aimed eventually at showing that requirement [2] of 
Lemma [5] (drift bounded from above) is satisfied on each locus = O. To this end we must first specify 
the various poles involved in the different Oi making up O = Oi A . . . A Ok and identify the region 
where = O (thus, "localizing"). For fixed £ > there can only be finitely many maximal linear segments 
^1, • • • , C 9D of length at least e. Fix a circle C centred in D of radius R, and locate the poles p^'", p^ 
at the intersection points on C of the line defined by li, sign ± chosen according to orientation (Figure [l]). 
Here R must be chosen to be large enough to fulfil asymptotics given below at ()2ip : in addition we require 
that R > diam(D) cscfcf)), where Scf) is the minimum modulo n of the non-zero angles between lines li, £j. 
The rationale for this is that if no other line segments are parallel to a given £1, and if x, y g £1, then the 
term corresponding to li {via the poles pf) in the minimum (fT9|) is then the unique minimizer. For example 
Vp+ 5 s {x,y) is the unique minimizer if k(x — ij, ^^-i-^ ^vt) is the unique minimum of the corresponding inner 
products. Hence R > diam(D) csc{4)) suffices to localize in the absence of parallelism, by a simple geometric 
argument indicated in Figure [21 This argument uses the remark (established by calculus) that 

exp (^-l||2i±v_p + ||^(^_y,!i±l_p+) 

is the minimum if no other 'pf is separated from D by the perpendicular to ii at p^ (which follows from the 
constraint of (fl6ll , which we have required for all poles p). 

In case some of the £i are parallel, then we must use the particular features of Vp_6_s as opposed to 
Vp_5, taking into account the expression (fl!]) for k. Suppose that x, ij G £i and li is parallel to £j. Establish 
cartesian coordinates (u, v) centred at ^^y^, for which £1 lies on the u-axis. Without loss of generality suppose 
that £j lies on the locus v = h. Suppose the circle C is centred at (uo,Vo) lying between £1 and Ij (Figure [3|. 

We need to show that, for all large enough R and all other poles pf for which (x — y, ^^y^ ~Pf ) 
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Figure 2: Let Scf) be the minimum modulo n of the non-zero angles between line segments £i, £j. If 
R > diam(D) csc(4)) and li is parallel to no other line segments then the poles , for a given li can be 
seen to supply the minimum in (fT9|) when x,y G Hi. 



same sign as — y 



' 2 



exp(-|ll^ 



exp(-l||^-pf||)(x-y,2S±v_p± 



> 



(20) 



and we now show that this will be the case if 1 < ff = S/R < 2 for large enough R (note that once R is large 
enough this is compatible with (fT6|) , which is our other requirement on R). 




Figure 3: Illustration of the movement of poles p^, when the centre of C is moved perpendicularly to 
the parallel line segments £i , Ij. 

First note that it suffice to consider the case when Vq is as large as possible. For if the centre of C is 
moved say from (uo,Vo) to (uo,h.] then the pole 'pf is moved further out on £j, while the pole p^^ is brought 
closer in on ii (see Figure [3|. Calculus shows that 

uexp ( — |- VuM-v2 j 



is increasing in u for < u < Sy — - — ^ (true for all feasible moves when cr = S/R > 1, since 
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R > diam(D) cscfcf)) and < <p < 7t/3), so it follows that such a move will increase the denominator and 
decrease the numerator of ([20|l . 

We now establish asymptotics which are uniform in (ucVq] e D; indeed, uniformly for — diam(D) ^ 
Uq < diam(D) and holding cr — S/R fixed, 

exp(-|||^-p,+ ||)(x-y,i±^-p+) ^ VR^^ + uo exp(-l(VR^^ + Uo)) 



exp(-l||^-pf||)(x-y,2S±v_p±) " R + Uo exp (-1 + ^oP + 

= " ^ u^ — ) ^ ^^P ~ - + t\/(R + uo)^ + - (R + uo)) 

^ (l - R^(^ + oi^-'))) X exp (^1 + o(R-^) + ^(.^^ + o({R + Uo)-^)); 

- + (21) 

It follows from these uniform asymptotics that if 1 < cr < 2 then ([20|) is satisfied for all large enough R. 

As a consequence of these considerations, and bearing in mind the continuity properties of the criterion 
ratio on the left-hand side of ([20|) , if we fix S/R € (1,2) and R large enough then (noting ||x — y\\ > e) it 
follows that for all sufficiently small ri > 0, if x, y e D n (£i © ball(0,ri)) and one of x, y G 9D then the active 
component of (fT9|) is the one involving li via the poles -pf. 

The final step is to establish Lemma [5] requirement [2] (drift bounded from above). When is the 
active pole, the drift is given by 

X-Y+k(X-p±) -^(X-Y, Z)e,v(X))dL'^ + ^Y-X - k(Y-p±)-^(X-Y, Z>e,T{Y))dL^ 

+ (n- |tr(JI + J^))dt 

(using the notation of the proof of Theorem [8}. In order to establish requirement [2] of Lemma [H it suffices 
to show non-positivity of 

(X - Y + K(X - p±) - J|(X - Y, Z) e,^(X) 

when X G 9D, and of 

Y - X - k(Y - p±) - ^(X - Y, Z> e,T{Y) 
when Y G 9D. Bearing in mind that ||X — Y|| ^ diam(D) and k ^ 6, if 

6 < X :r-^ (221 

diam{D)+R l + ^diam(D) 

then this follows directly when (X — Y,-v{X)) > £, (in case X G 9D) or when {Y—X,^/[Y]) > £, (in case Y G 9D). 

It remains to consider the case when this does not happen. By choosing f, small enough, we may then 
ensure that the two Brownian motions are both close to the same segment portion of the boundary. For the 
function 
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is continuous on the closed and bounded (and therefore compact) Euclidean subset 

H — {(x,'y,-v) : x e D,-li e 9D, ||x — yH ^ £,"v normal at y} , 
and vanishes only on Uiil^iyi"^) G H : x,y e Consequently for any r] > we can find £, > such that 

(y-x,-v(y)) < ||x-y||f, 
forces X, i| e © ball(0,ri) for some i. If r| is chosen as above then this in turn forces to be the pole for 

We now argue for the case when y G 9(D); the case when x e 9(D) is similar. We can choose r| > as 
small as we wish: we therefore require that ^l-TiV£2 -Ti/e > 0. 




n 



Figure 4: Geometric construction underlying analysis of singular drift when y G 9D and (y — x,'v(y)) < 

Suppose then that ij is at perpendicular distance h, < t] from the line through li (FigureH]). If x is further 
away from Ix than y, then the singular drift will certainly be non-negative once ([22]) is satisfied, since the 
segment y — x then makes a smaller angle with y[vi) than does either y — pf or ^^-^-^ • Otherwise (using 
the angle notation indicated in Figure [4]) we require non-negativity of 

£sin(a— |3) — 6 X (diam(D) + R)(l + ^ diam(D)) sin a, . 

Hence we require 

^ £ sin(a— 13) 

(diam{D) +R)(1 ^diam(D) sina 

£ 

— -. (cos 6 — sin 6 cot a) 

(diam(D) + R)(l + ^diam{D))^ ' ' ^ 

^ £ 

^ (diam(D) R)(l + ^ diam(D)) 

^ £ 

^ (diam(D) R)(l ^ diam(D)) 

Thus for all sufficiently small 6 and sufficiently large R the singular drift is non-positive (the case of x e 9D 
following by the same arguments) and so the drift is bounded above as required, thus completing the proof. 

□ 




U^-i) • (23) 
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4 Conclusion 



The above shows how rather direct potential-theoretic methods permit the deduction of non-existence of 
shy couplings for reflected Brownian motion in bounded convex domains, so long as either (a) the domain 
boundary contains no line segments (Theorem!?} or (b) the domain is planar (Theorem[8]). One is immediately 
led to the conjecture that there are no shy couplings for reflected Brownian motion in any bounded convex 
domains. In the planar case one can extend the substantial linear portions of the boundary to produce a 
finite set of poles p^, leading in turn to the key function CD expressed as a minimum of a finite number 
of simpler functions. Similar constructions will in the general case lead to a continuum of possible poles, 
contained in an essential continuum of different hyperplanes, and this variety of poles is an obstruction to 
generalization of the crucial localization analysis in the proof of Theorem [H Careful geometric arguments 
allow progress to be made in the case of bounded convex polytopes, for which attention can be confined to 
a finite number of hyperplanes and hence to sets of poles forming hypercircles via intersection of a large 
hypersphere with the hyperplanes. However even in this limited special case tedious arguments are required 
in order to overcome problems arising from intersections of the hypercircles, and we omit the details. 

Further variants on the general theme of s hy-ne ss are possible. For example, it is not hard to use the 



comparison techniques described in ljost et al. ( 19981 Chapter 2) to show the following: if M is a Riemannian 



manifold with sectional curvatures all bounded above by a finite constant k^, then there exist £-shy couplings 
of Brownian motion on M for all £ satisfying 

£ < min{^, injectivity radius of M} . 

These couplings are geometric versions of the perverse coupling described in Sec tion [TJ 



However the major challenge remains the conjecture of Beniamini et al.l (j2007l . Open problem 4.5(ii)), who 



ask whether there is any simply connected planar domain which supports a shy coupling of reflected Brownian 
motions. The present work br i ngs us close to understanding shy coupling for convex domains; progress in 



resolving the iBeniamini et al.l (|2007l ) conjecture requires development of completely new techniques not 



dependent on convexity at all. 
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